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ABSTRACT 



The A-type star HR 6412 = V2368 Oph was used by several investigators as a photometric comparison star for the known eclipsing 
binary U Oph but was found to be variable by three independent groups, including us. By analysing series of new spectral and 
photometric observations and a critical compilation of available radial velocities, we were able to find the correct period of light and 
radial-velocity variations and demonstrate that the object is an eclipsing and double-lined spectroscopic binary moving in a highly 
eccentric orbit. We derived a linear ephemeris r,^i„.i =HJD (2454294.67±0.01) -I- (38^327 12±0'?00004) xE and estimated preliminary 
basic physical properties of the binary. The dereddened UBV magnitudes and effective temperatures of the primary and secondary, 
based on our light- and velocity-curve solutions, led to distance estimates that agree with the Hipparcos distance within the errors. 
We find that the mass ratio must be close to one, but the limited number and wavelength range of our current spectra does not allow a 
truly precise determination of the binary masses. Nevertheless, our results show convincingly that both binary components are evolved 
away from the main sequence, which makes this system astrophysically very important. There are only a few similarly evolved A-type 
stars among known eclipsing binaries. Future systematic observations and careful analyses can provide very stringent tests for the 
stellar evolutionary theory. 

Ksy words, stars: early-type - stars: binaries - stars: A - stars: individual: V2368 Oph , U Oph 



1. Introduction 

HR 6412 = HD 156208 has often been used as the photometric 
comparison star for a we ll-known eclipsing binary U Oph, which 
exhibits apsidal motion (Huffer& Kopal 1951; Koch & Koeglet 
il977; Wolf et al. 2002; Vaz et al. 2007). McAlister et al. ( 1987) 
reported that HR 6412 is a speckle-interferometric binary with a 
separatio n of 0'.'136 and es timated the orbital period to 72 years. 
However. iMcAlister et al.l ([T993) could not resolve this pair and 
concluded that the original detection had been spurious. 

In Table [1] we summarize the various determinations of the 
yellow magnitude of V2368 Oph published by several authors. It 
seems to indicate that no secular variations in its brightness have 
been recorded, since the scatter of values of a few hundredths of 
a magnitude is quite normal for yellow magnitudes recorded in 
different photometric systems. 



Table 1. Published yellow magnitudes of V2368 Oph. 



Mag. Source 



Photometric system 



6"? 17 
6"? 19 
6"? 16 
6"?178 
6"? 177 
6"?18 



Easen (1955) 
Stokes (1972) 
Becker etal. (1975) 
Gr0nbech & Olsen (1976 ) 
Sowell & Wilson (1993) 
van Gent (1982) 



(P, V)e system 

uvby and Hf3 

Cousins' values 

uvby; 1965 - 1970 

uvby; Nov. 1988 & Apr. 1991 

BVR; 1970 



Send ojfprint requests to: Petr Harmanec 
e-mail: : PetrHarmanec@mff.cuni.cz 

* Based on new spectral and photometric observations from the 
following observatories: Dominion Astrophysical Observatory, Hvar, 
Ondfejov, San Pedro Martir, Tubitak National Observatory, and ASAS 
service 

** Tables 2, 3, and 4 are available only in electronic form at the 
CDS via anonymous ftp to cdarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qc at ?J/A-hA/| 



During our 2001 observations of U Oph at San Pedro Martir 
Observatory (SPM hereafter), HR 6412 was also used as the 
comparison star We noticed large changes in this comparison on 
JD 245207 1 .7 1-1 .85. Upon a literature search , we found that the 
variab ility of HR 6412 has been discovered byl Perrvman & ESAl 
(Il997h . who classified it as an eclipsing binary with a period 
of 7'?70. iKazarovets et al.l (1 19991) then assigned it the variable- 
star name_V2368_Oph. The variability has also been confirmed 
by I Vaz et al. (200 7|), w ho mention that the period found by 
iPerrvman & ESA (119971) was incorrect but give no other details . 
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Fig. 1. One complete Ondfejov red spectram of V2368 Oph, 
taken on HJD 2454357.2917, which shows that the only stronger 
lines, suitable for the RV measurements, are Ha and the 
Sill 6347 and 6371 A. The spectrum contains many water 
vapour and oxygen telluric lines and the interstellar Une at 
6613.56 A. 



The main goal of this study is to publish the first correct and 
accurate linear ephemeris of V2368 Oph, which can be used to 
correct earlier photometric observations of U Oph. We also de- 
rive preliminary orbital and light-curve solutions and show that 
they can lead to self-consistent basic physical properties of the 
binary. However, considering the limited amount and hetero- 
geneity of our observational material, we do not aim to deter- 
mine the final, accurate physical elements of the system. 



2. Observational material used 

2.1. Photometry 

After realising that V2368 Oph is a variable star, we started sys- 
tematic UBV observations of it at Hvar and SPM observatories. 
A limited set of UBV observations was also obtained by HA 
at the Turkish National Observatory. Besides, we compiled the 
Hipparcos H„ obse rvatio ns and V photometry from the ASAS 
project (Poimanskill2002h . Details on data sets and their reduc- 
tion are in Appendix lAl and all individual UBV and V observa- 
tions with their HJDs are provided in Table 2 (electronic only). 



2.2. Spectroscopy 

Simultaneously with photometric observations, we also begun 
to collect electronic spectra in Ondfejov, San Pedro Martir, and 
Dominion Astrophysical Observatory (OND, SPM, and DAO 
hereafter). A detailed discussion of all spectra, their reduction, 
and a journal of observations can be found in Appendix iB] 

Here, we only want to add a few comments relevant to fur- 
ther analyses. The only spectral region that is available in the 
spectra from all three observatories is the red region containing 
only three strong enough spectral lines suitable for the RV mea- 
surements: the Balmer Ha line and the doublet of Si II 6347 and 
6371 A lines (see Figs. [T] and |2]i. The SPM spectra also cover 
the region of Mg II 448 1 A line, in which both components are 
clearly seen, so this line was also found suitable for the RV mea- 
surements. The RV measurements were carried out in three dif- 
ferent ways. (1) We used the program SPEFO dHorn et al.i[l996t 
ISkodal[r996l) . which permits the RV measurements via sliding 
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Fig. 2. The Ha and Sill 6347 and 6371 A fine profiles from 
the Ondfejov spectrum taken near one periastron passage on 
HJD 2454636.4608. One can clearly see the lines of both binary 
components, separated by more than 200 km s"' in their radial 
velocity. 



direct and flipped line profiles on the computer screen until a 
perfect match is obtained. (2) We fitted the observed line pro- 
files by two Gaussians shifted in their positions in such a way 
as to obtain the best match of the observed, often blended line 
profiles of the primary and secondary. (3) We verified our RV 
values by an automated fitting of a combination of two synthetic 
spectra — selected from the Ondfejov library of synthetic spec- 
tra, prepared and freely distributed by Dr J. Kubat — to each of 
th e observed spectra using a simplex algorithm. The (defined 
in iBroz et al ]f2010) was calculated for the entire red spectrum 
in the wavelength range 6256 to 6768 A, which includes all the 
individual lines analysed previously. 

The second method is preferred as it returns the most Ukely 
velocity amplitudes. The results of the third method are statisti- 
cally compatible. The first method (SPEFO RV measurements) 
usually underestimates the true semi-amplitudes of the orbital 
motion, because of the line blending, especially for the steep 
wings of Ha. We, therefore, used the SPEFO Ha RVs only for 
the initial search of the orbital period, to combine them with 
older published RVs, measured in a standard way from the pho- 
tographic spectra. 

For a ll red spectra, we followed the procedure outlined in 
iHom et al. ( 1996) and measured RVs of a selection of unblended 
telluric lines in SPEFO. We then used the difference between the 
calculated heliocentric RV correction and the true mean RV of 



2 



p. Harmanec et al.: Eclipsing binary V2368 Oph 



> 

OS 




0.4 0.6 
orbital phase 

Fig. 3. The radial-velocity curve of the primary component of 
V2368 Oph based on He RVs measured in SPEFO for our OND, 
DAO, and SPM spectra (filled circles) and published RVs (open 
circles) for the period of 38'?3307 from the FOTEL orbital solu- 
tion with phase zero at periastron passage. See the text for de- 
tails. 

telluric lines to correct the zero point of the RV scale individ- 
ually for each spectrogram. Regrettably, these corrections were 
less reliable for the majority of the red SPM spectra, which only 
contain very weak telluric lines owing to the high altitude of that 
observatory. No such corrections were possible for the blue SPM 
spectra in the absence of telluric lines, of course. 

Individual Ha, Si II, and Mg II RVs, measured in a standard 
way in SPEFO and via Gaussian fits to line profiles with their 
HJDs of mid-exposures are in (electronic only) Table 3. The rec- 
tified and wavelength-calibrated spectra are in (also electronic) 
Table 4. 



3. Preliminary analysis and search for a correct 
orbital period 

As our observations progressed, it soon became obvious that the 
or bital period must be much longer than the 7'?70 period reported 
by|P errvman &ESAl(ll997h and that the orbit had to have a high 
eccentricity, since we were observing a constant brightness and 
only small RV changes. When we finally succeeded in observ- 
ing a decline into the minimum on the night JD 2454294.35- 
4.54, we were able to combine it with earlier minima recorded 
by Hipparcos, ASAS, and our discovery observation at SPM and 
to obtain the first guess that the period should be close to 38 
days. Continuing spectroscopic observations then allowed us to 
cover parts of two periastron passages on JD 2454366-67, and 
JD 2454636 and an iterative analyses of the RV and light curves 
allowed us to estimate the value of the orbital period to 38'?33. 

Ther e are t wo limited sets of earlier RV measurements. 
IChristid (119251) pubUsh ed five RVs covering the interval 
JD 2423233.8-995.8, and lPalmer et"al](ll968l) published another 
five RVs from low-dispersion spectra covering JD 2437441.6- 
740.7. We combined these RV s with our own RV measure- 
ments in the program SPEFO (iHorn et al.lll99 6; Skoda" 1996t) 
for the Ha line and used the FOTEL program (Had raya 1990l 
l2004ah to derive preliminary orbital elements and a more accu- 
rate value of the period. We obtained P = 38^3302 + ^0015, 
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Fig. 4. The final radial-velocity curves of both binary compo- 
nents of V2368 Oph based on Gaussian fits (circles) and the the- 
oretical RV curves based on our PHOEBE solution (lines). The 
lower panel shows the 0-C residua from the solution. Typical 
observational uncertainties are 10 km s Orbital phases from 
ephemeris(|2|l are used. 



model fit to the data per 1 observation being 10.3, and 4. 1 km s"' 
for the old and new RVs. 

Since the narrow and steep photometric eclipses are very 
sensitive to the phase shifts, we used an interactive computer 
program (written by HB), which allows the user to display the 
phase diagrams based on the observed data in the neighbourhood 
of the eclipses for various smoothly varied values of the orbital 
period. This way we found that the true orbital period must be 
very close to the value of 38^3272. 



4. Towards basic physical properties of the 
binary 

To obtain self-consistent physical properties of the components 
and the binary system, we had to proceed in an iterative way. We 
selected several stronger lines seen in both binary components 
and derived their RVs via Gaussian fits to line profiles. In par- 
ticular, we used the Sill 6347&6371 A and Ha lines, available 
in all spectra, and Mg II 4481 A, measurable in the SPM spec- 
tra. For Ha, the Gaussian profiles were not optimal so we ten- 
tatively disentangled t he Ha profiles, using the KOREL program 
(lHadravalll995l Il997l |2004b. 2051 and used the disentangled 
profiles instead of GaussiansQ 

We alternativel y used the p rograms PHOEBE 

(|Prsa&Zwittedl2005l l2006l) based on the IWilson & DevinnevI 



e = 0.551 + 0.010, w = 355.2 + 1.6, Ki = 59.71 ± 4.1 km s"'. 
Told - 17.0 ± 3.4, and ynew - 10.07 + 0.50, the rms errors of the 



' It would seem logical to derive the orbital solution directly with 
KOREL. However, due to heterogeneity of the available spectra, their 
different spectral resolutions and relatively limited number, this proce- 
dure was not satisfactory in the given case. 
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Table 5. Published uvby and li/3 observations of V2368 Oph. 



V 


b-y 






^uvby 


H/3 




Source 




0'^167 


o";o90 


i";i95 


5 


2.850 


4 


Crawford et al. (1972) 


6'-^ 19 


0"?177 


0'';095 


1™216 


3 


2.860 


3 


Stokes (1972) 


6™178(1) 


0"?161(2) 


0";095(2) 


r;204(2) 


2 






Gr0nbech & Olsen (1976) 












2.870(2) 


3 


Gronbech et al. (1977) 


6";i77(5) 


0"? 189(2) 


0";072(5) 


1™203(5) 


3 






Sowell& Wilson (1993) 



Table 6. The final combined light-curve and RV-curve solutions obtained with PHOEBE. 



Element 




Primary 


Binary 


Secondary 


Primary 


Binary 


Secondary 


P 


(d) 




38.327115(43) 






38.327118(43) 




^peiiaslr. 


(RJD) 




54291.039(11) 






54291.042(11) 




^min.I 


(RJD) 




54294.670(17) 






54294.668(17) 




^min.II 


(RJD) 




54287.498(17) 






54287.495(17) 




e 






0.51527(14) 






0.51524(14) 




CO 


n 




359.33(20) 
86 165(^??1 






359.41(20) 




r 




0.0473 




0.0460 


0.0477 




0.0464 


n 




23.29(11) 




24.75(12) 


23.05(11) 




23.71(11) 


a 


(Ro) 




83.67(67) 






83.70(68) 






(kms-') 


65.7(8) 




62.9(9) 


64.3(6) 




64.3(6) 








1.044(15) 






1.0 fixed 




Tsff 


(K) 


9300 fixed 




9500(200) 


9300 fixed 




9500(200) 


M 


(Mq) 


2.62(2) 




2.74(7) 


2.68(8) 




2.68(8) 


R 


(Ro) 


3.96(2) 




3.84(2) 


3.99(3) 




3.87(2) 


Mboi 


(mag) 


-0.31(9) 




-0.34(9) 


-0.33(9) 




-0.36(9) 


logg 


[cgs] 


3.66(1) 




3.71(1) 


3.66(3) 




3.69(3) 




l/band 


0.5029(30) 




0.4971 


0.5024(30) 




0.4976 




fiband 


0.5004(32) 




0.4996 


0.4995(32) 




0.5005 




f/band 


0.4942(37) 




0.5058 


0.4927(38) 




0.5073 


V 


(mag.) 


6.913(16) 


6.167(11) 


6.926(16) 


6.913(16) 


6.166(11) 


6.924(16) 


B 


(mag.) 


7.137(17) 


6.386(11) 


7.139(17) 


7.139(17) 


6.385(11) 


7.137(17) 


U 


(mag.) 


7.363(20) 


6.598(13) 


7.338(20) 


7.367(20) 


6.599(13) 


7.335(20) 


Vo 


(mag.) 


6.26(8) 




6.29(8) 


6.26(8) 




6.29(8) 


(B-V)o 


(mag.) 


0.021(11) 




0.014(11) 


0.021(11) 




0.014(11) 


(U-B)o 


(mag.) 


0.077(15) 




0.054(15) 


0.078(15) 




0.053(15) 


No. of obs. 


UBV/RW 




1913/268 






1913 /268 










1726 






1730 





Notes. Columns 2^ contain a solution based on the free convergence; columns 5-7 contain the solution for the fixed mass ratio of 1 . All epochs are 
in RJD = HJD - 2400000. Probable elements and their formal error estimate s are provided, where CI is the value of the Roche-model potential used 
in the WD program, and Lj (j = 1,2) are the relative luminosities of the components in individual photometric passbands. They are normalized so 
that Li + L2 = 1. The number of observations represents the sum of RVs of the primary and secondary and a sum of individual observations in all 
passbands. 



(Il971h program, and FOTEL, akeady used in the first step, 
to derive preliminary values of some critical parameters. To 
obtain the best possible estimate of the RV semi-amplitudes, 
we allowed for individual systemic velocities for each of 
the ions used. In particular, we found the systemic velocities 
of 2.97+0.55, 0.14+0.54, and 9.9+6.0 for Sill, Mgll, and 
Ho-, respectively. The differences between these values are 
probably insignificant considering their errors and the inability 
to check the zero point of the RV scale for the blue spectra via 
measurements of the telluric lines. Since PHOEBE can treat 
only one joint systemic velocity, we subtracted the values of 
respective systemic velocities from the observed RVs and used 
these shifted RVs from all three ions as one dataset for the 



primary and another one for the secondary in PHOEBE. We then 
naturally kept the systemic velocity fixed at zero in PHOEBE 
solutions. 

In the latest (development) version of PHOEBE that we are 
using, the convergence is governed by minimization of a cost 
function defined in the case of our datasets as 



X 



1 



(1) 



;=I 



where index p denotes the individual photometric passbands, 
CTp their standard deviations per 1 observation, Np is the num- 
ber of individual observations for p-\\\ passband, w, are standard 
weights of individual observations, and /j and Sj the observed 
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Fig. 5. The observed V-band light curve compared to the theo- 
retical one, based on our PHOEBE solution. The lower panels 
show a zoom of the curves in the neighbourhood of both bi- 
nary eclipses and the 0-C residua from the model fit. Typical 
1-cr observational uncertainties are O'^l. Orbital phases from 
ephemeris (|2|i are used. 



the binary between 8900 K, and 9400 K and a mean log g be- 
tween 3.50 and 3.61. 

Since both stars are detached well even near periastron, the 
light-curve solution basically does not depend on the exact value 
of the mass ratio. Considering this, a preliminary PHOEBE solu- 
tion was used to derive relative luminosities of both components, 
UBV magnitudes of the binary at maximum light from the UBV 
observations transformed to the standard system and from them 
the UBV magnitudes of the primary and secondary in each pass- 
band. These were then dereddened in a standard way, assum- 
ing Ay - 3.2 E(B-V). The dereddened magnitudes and indices 
confirmed the spectral type of A2 for the primary. Using Flower| 
(1996) calibration of (B-V)q indices vs. bolometric corrections 
and TefF, we estimated the effective temperature of the primary 
to (9300+200) K. Keeping the value of 9300 K fixed, we de- 
rived a freely converged PHOEBE solution, which is presented 
in detail in Table |6] The rms errors of all converged parame- 
ters in PHOEBE are derived from a covariance matrix. For other, 
deduced parameters, we propagated the errors to obtain the esti- 
mates given in the Table. Since a realistic error of the effective 
temperature of the primary is about +200 K, this must imply that 
the formal error of the effective temperature of the secondary, 
estimated from a covariance matrix in PHOEBE, is too low and 
must also be about ±200 K. 

It is clear at first sight that this solution is not satisfactory 
since it leads to a model in which the more massive component 
is the less evolved of the two. We believe that the problem lies 
in the limited quality of our radial velocities. As a matter of fact, 
the solutions for RVs of individual ions oscillated between 0.95 
and 1.05 in the resulting mass ratio. At the same time, since 
both binary components are well detached, the light curve so- 
lution is stable and basically does not depend on the mass ratio. 
Therefore, we derived another PHOEBE solution, this time with 
the mass ratio fixed to 1.0. This solution is also provided in the 
last three columns of Tableland we take it as the reference solu- 
tion for the following discussion. The corresponding RV curves 
and the light curve in the V band are shown in Figs. ID and|5] 
respectively. It is seen that there is little difference in all parame- 
ters, which do not depend on the mass ratio between the two so- 
lutions. For completeness, we also derived another solution for 
the mass ratio of 0.95. This led to a slightly worse - 1770, 
but the photometric elements were again very similar to those 
two shown in Table |6] 

The solutions led to the following linear ephemeris, which 
should enable correction of existing photometric observations of 
U Oph secured differentially relative to V2368 Oph: 



and calculated fluxes, respectively. The value of the function 
is tabulated along with the solutions. 

Although it should be possible to derive the effective tem- 
peratures of_botii_binarycorin)on^^ calibrated UBV pho- 
tometry ('Prsa & Zwitter"2006VWilson 2008), the propagation of 
errors often leads to unreliable results. For that reason we re- 
stricted ourselves to the standard approach of estimating the ef- 
fective temperature of the primary from the dereddened colours 
and observed spectra, then keeping its value fixed in the solu- 
tions. 

There are four sets of uvby observations and three sets 
of H/3 photometry of V2368 Oph - see Table |5] Using 
the program UVBYBETA written by TT Moon and modified 
by R. Napiwotzki, whi c h is based on a calibration devised 
by iMoon & Dworetskvl (Il985l) . we found that those sets of 
Stromgren photometry imply a mean effective temperature of 



r„i„.i = HJD 2454294.67 + 38^327 12 xE . (2) 

It is encouraging to note that a separate dereddening of the 
UBV magnitudes of the primary and secondary led invariably 
to E(B-V) = 0"!20 and to distance moduH of 6.46(12), and 
6.47(12) for the primary and secondary, respectively. The dered- 
dened values of the secondary indicate a slightly earlier spec- 
tral type, in accordance with its higher effective temperature ob- 
tained from the PHOEBE solution. We note that E{b-y) derived 
with the program UVBYBETA from the uvby values of Table|5]is 
0™15, which agrees well with the E(B-V) derived by usQ 

The parallax of V2368 Oph was obtained by the 
ESA Hipparcos mission , and its originally published value 
(iPerrvman & ESA 1997) is 0'.'00554 +^.'00086, while an im- 
proved value obtained by van Leeuwen (.2007 a^h) reads as 

- Note that E(b - y) = 0.74£(B- V). 
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(X' 00455 +a.' 00048. The distance modulus obtained from our 
photometric solution implies a parallax of 0'.'00506, in excellent 
agreement with the above values, deduced from the Hipparcos 
observations. 

5. Stellar evolution of the components 

Since V2368 Oph is a detached binary with no evidence of 
mass transfer, one can use a one-dimensional program for stellar 
evolution to see whether the observed properties of the binary 
agree with the model prediction. To this end, we used the stellar- 
evolution module MESAstar bv lPaxton et all (1201 ih . 

We first calculated the model evolution for the masses Mi - 
2.62 Mo and M2 = 2.74 Mq which follow from the free PHOEBE 
solution with the lowest for the combined photometric and 
RV data (Table |6] left). We assumed the same metallicities of 
Z = 0.02 for both components, of course, the helium abundance 
Y = 0.28 and the mixing-length parameter a = 2.0. The result is 
compared to the observed binary properties in the Hertzspmng- 
Russell diagram, and the T^ff vs. radius diagram (Figure |6] left 
panels). Even though there are uncertainties in temperatures, lu- 
minosities and masses of the individual components (refer to 
Table|6l), their differences are established much more accurately; 
e.g., the difference in temperatures T2 - Ti ^ 200 K is always 
present in PHOEBE solutions since this is enforced by the ob- 
served light curve and colour indices. 

There is a clear disagreement between the photomet- 
ric/spectroscopic observations and the predicted stellar evolution 
in this case. Since the mass ratio q = M1/M2 - 0.96 differs 
significantly from 1, the calculated luminosities of the compo- 
nents are always very different owing to a strong dependence 
of the stellar evolution on the mass (log L2/LQ reaches ^2.0 and 
logLi/Lo ^ 1 .9), while the observed luminosities are rather sim- 
ilar (logLi/Lo ^ \QgL2ILQ ^ 2.03). A change in neither metal- 
licity nor in the mixing-length parameter could alter this result 
since a different value of Z would shift both tracks in the same 
direction, and varying a from 1 .5 to 2.5 does not alter evolution- 
ary tracks significantly before the red giant branch is reached. 

As a second test, we took the mass ratio q close to 1, which 
is still compatible with the photometric/spectroscopic observa- 
tions from a statistical point of view (Table |6] right). Because 
the stellar evolution is very sensitive to the stellar mass, we 
may actually use this approach to constrain the mass ratio of 
V2368 Oph. The most sensitive indicator seems to be the tem- 
perature - there is approximately a 200 K difference between 
Ji and T2, which corresponds to a 0.004 to 0.008 Mq differ- 
ence between M\ and M2, according to our tests. If we use 
Ml = 2.682 Mq and M2 = 2.678 Mq we also have to decrease 
the metallicity to Z = 0.014, which shifts both the evolutionary 
tracks towards higher T and L, in order to match the observed 
state of V2368 Oph (Figure|6] right panels). Another possibility 
would be to slightly increase the masses to Mi = 2.760 M© and 
M2 - 2.752 Mq and to retain the Z = 0.02 value. To conclude, 
it is possible to find a consistent solution for the available pho- 
tometry and spectroscopy and the stellar evolution, even though 
the parameters presented above cannot be considered as final, 
because the total mass of the system is not yet constrained pre- 
cisely enough. 

From the standpoint of stellar evolution, V2368 Oph is a 
very interesting evolved system with both components leaving 
the main sequence. It is in a rapid phase of evolution and conse- 
quently may serve as a very sensitive test case for the stellar- 
evolution programmes, provided new, accurate RVs and pho- 
tometric observations are acquired. Considering the relatively 



short distance of the binary from us, its interferometry would 
also be of utmost importance, providing the angular separa- 
tion of the components and orbital inclination, consequently a 
much more precise parallax and the total mass of the system. 
Interferometry can also provide independent constraints on the 
component radii. 

6. A comparison with synthetic spectra and the 
rotation of the binary components 

As another consistency check, we disentangled t he Si II and Hg 
line profiles with th e help of the KOREL program (iHadraval 19951 
Il997 ','2004b','2005'). keeping the orbital parameters from the sec- 
ond PHOEBE solution fixed, but using the mass ratio of 0.998 
(considering the discussion above). In Fig. |7] the disentangled 
line profiles, normalized to their individual continua using the 
relative luminosities derived by PHOEBE, are compared with 
the synthetic line profiles from the Ondfejov library of synthetic 
spectra prepared^and freely distributed by Dr J. Kubat — see, 
e.g. iHarmanec et al.l (fl997ai) for details. We used the synthetic 
spectra for the parameters close to the PHOEBE results, namely 
Teff = 9500 K and log g = 3.5, rotationally broadened in SPEFO 
to 40km s"' for the primary and to 90 km s ' for the secondary. 
Varying these values for more than +5 km s"' would result in 
a significant disagreement in the depths and widths between the 
observed and synthetic line profiles. 

The agreement between the observed and synthetic spectra 
is satisfactory, lending some credibility to our result. We warn, 
however, that the heterogeneity of our material means that there 
can still be rather large uncertainties in the derived masses, radii, 
and luminosities. Another study based on rich and homogeneous 
observational material is therefore desirable. 

Taken at face value, the observed projected rotational veloc- 
ities and the PHOEBE solution would imply the stellar rotational 
periods of 5^0 and 2^2 for the primary and secondary, respec- 
tively, while the spin-orbit synchronization in periastron would 
imply a rotational period of 10^5. It is notable that in their de- 
tailed study of ano ther evolved A-type binary with eccentric or- 
bit, 6^ Tau, Torres et al.l (1201 ll) also found the projected rota- 
tional velocity for the secondary roughly twice as high as for the 
primary. 

In passing, we wish to mention that it is also possible to 
make a theoretical prediction of the interna l structure co nstant 
based on the current evolutionary models of IClaretl (12004;) with 
the standard chemical composition of iX,Z) = (0.70,0.02): 
logA:2 = -2.498 for both binary components. Taking the val- 
ues of the eccentricity and fractional radii from Table |6]into ac- 
count, we can predict a very slow apsidal-motion rate of Wobs = 
0.00028 deg/cycle, which is only 0.27 deg/century. The rela- 
tivistic contribution to the apsidal motion is substantial: (ii,ei = 
0.00020 deg/cycle or about 70% of the total apsidal-advance rate 
(Gimenez 1985). In other words, there is little chance of detect- 
ing a measurable apsidal motion for this binary in the foreseeable 
future. 
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Fig. 6. Left panels: A Hertzsprung-Russel diagram (top), and the Teff vs. radius diagram (bottom) showing the stellar evolution of the 
primary and secondary components of V2368 Oph. It was computed with the MESAstar module, for the masses Mi - 2.62 Mq, 
M2 - 2.74 Mq and for the metallicity Z = 0.02. The evolutionary tracks are plotted by thick lines from ZAMS up to the age 
3.877 X 10** y. The ranges in temperatures T[, T2, luminosities Li, L2 and radii Ri, R2 inferred from photometry /spectroscopy are 
denoted by lines (refer to Table |6]l. There is a strong disagreement between the observations and the stellar evolution, especially 
for the primary. Right panel: same for the masses Mi = 2.682 Mq, M2 = 2.678 Mq (i.e., the mass ratio very close to one), and a 
different value of the metallicity Z = 0.014. The thick lines are terminated at the age 3.938 x 10*^ y. 
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Fig. 7. A comparison of disentangled (solid lines) and synthetic (dotted lines) profiles of the Si II fine (left panels) and Ha Une 
(right) for the primary (top panels) and secondary components (bottom). See the text for details. 
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Appendix A: Photometry 

Here we provide detailed comments on the photometric obser- 
vations used in this study and the way we treated them. 

1. Hvar observations were secured in 2007, 2008, and 2009 
and reduced and transformed into the standard UBV sys- 
tem with the HEC22 release 16.1 reduction program via 
nonlinear tran sformation formulae (iHarmanec et al.l Il994i 
iHarmanec & H orn 1998). This recent version of the program 
allows modelling of variable extinction during the observing 
night, which significantly improves the accuracy of the ob- 
servations. The typical rms errors of the multinight fit to all 
standard stars used to define the transformation formulae in a 
given observating season are 0"X)08 for V and B, and 0"010 
for U. This is similar for observations from other observing 
stations reduced with the help of HEC22. 

2. San Pedro Martir observations were collected during the 
observational runs in 2001-2003 and 2007. Observations 
were reduced and transformed into the standard UBV sys- 
tem with the HEC22 release 14.1 reduction progra m via 
nonlinear transformation formulae ( Harmanec et al.l Il994t 
IHarmanec & HornllT998l) . At that station V2368 Oph was 
used in 2001 as a recommended comparison star for the 
eclipsing binary U Oph and its magnitude difference rela- 
tive to several comparison stars (HD 183324, HD 187458, 
HD 161132, HD 153808, and HD 144206) was derived. 
For all these stars, save HD 183324, the magnitudes and 
colours are well established from the calibrated Hvar all- 
sky photometry. HD 183324=V1431 Aql was found to b e 
a small-amplitude A Bootis variable ( Kuschnig et alj[l"994l) . 
It served for a long time as a comparison star for observa- 
tions of V923 Aql and VI 294 Aql in the Photometry of the 
Bright Northern Be Star Programme (Tiarmanec et al.l ll982[ 
[1994; Pavlovski et al. 1997; Harmanec et al. 1997b) and its 
variability on longer time scales is safely excluded by nu- 
merous Hvar observations. The mean all-sky UBV magni- 
tudes of HD 183324 are accurately derived. It was actu- 
ally used as a comparison for V2368 Oph only on the night 
JD 2452065, and we feel that its 2 mmag microvariability 
is not critical for the purpose of this study. For 13 observa- 
tions secured on JD 2452071 (when we recorded the first 
eclipse of V2368 Oph), it was not possible to derive the dif- 
ferential values for it so we adopted its all-sky values in- 
stead, since enough standard stars had been observed during 
the night, and the nightly transformation coefficients (extinc- 
tion and its variations and the zero points) could be derived. 
As soon as we realized that V2368 Oph is a variable, its 
subsequent observations in 2002 and 2003 were carried out 
differentially, relative to HD 154660=HR 6361. This A9V 
star is a visual binary ADS 10347 A with a close compan- 
ion ADS 10347B at 20'.'3, which is some 3™35 fainter than 
ADS 10347 A. The 2007 observations were obtained with 
a larger diaphragm so that the light of the visual compo- 
nent ADS 10347B was recorded with the brighter compo- 
nent ADS 10347A = HD 6361. We carried out dedicated 
observations at Hvar to derive the total magnitude of both 
visual components and added this value to the magnitude 
differences var - comp. from this season. In all other in- 
stances, observations were carried out in such a way as to 
keep ADS 10347B outside the diaphragm. 

3. Hipparcos all-sky Hp broad-band magnitudes secured be- 
tween 1989 and 1993 (Perryman & ESA ( 1997)) were ti-ans- 
formed to the standard Johnson V magni tude with the non - 
linear transformation formula derived bv iHarmanecl (Il998l) . 



Table A.2. Comparisons and check stars of V2368 Oph 



HD/BD 


Other ident. 


V 


B-V 


U-B 


1 j'+uOU 


rlK DjOi 


D.J J / 


noil 


n 1 


-01° 3292B 


ADS 10347B 


9.71 


0.66 


0.15 




ADS 10347AB 


6.308 


0.227 


0.094 


154895 


HR 6367 


6.058 


0.075 


0.028 


183324 


V1431 Aql 


5.801 


0.086 


0.067 


187458 


HR 7550 


6.660 


0.426 


-0.056 


162132 


HR 6641 


6.494 


0.085 


0.075 


153808 


6 Her 


3.916 


-0.024 


-0.088 


144206 


V Her 


4.739 


-0.096 


-0.326 



Magnitude and colours of ADS 10347 AB are values resulting from co-added flux of HR 
6361 and BD -01° 3292B measured simultaneously though a larger diaphragm in the 
photometer. 

Table B.l. Journal of spectroscopic data of V2368 Oph 



Spg- 


Time interval 


No. of 


Station, telescope 


no. 


(HJD-2400000) 


spectra 


and instrument 


A 


54266.4-54638.4 


11 


OND 2.0-m, grating spg. 


B 


54339.7-55102.6 


19 


DAO 1.22-m, grating spg. 


C 


52745.8-52749.0 


21 


SPM 2.1-m, echelle spg. 


D 


54193.9-54200.0 


21 


SPM 2.1-m, echelle spg. 


E 


54250.8-54253.0 


12 


SPM 2.1-m, echelle spg. 



The rms error of the fit per 1 observations is 0"X)067. For 
the solutions, the transmission and the limb darkening coef- 
ficients for the //p passband were considered, however All 
data with error flags larger than 1 and one deviating point at 
HJD 2448661.4682 were omitted. 

4. TNO (Tubitak National Observatory) observations were se- 
cured during two consecutive nights in 2003 and were 
reduced and transformed into the standard UBV sys- 
tem with the HEC22 release 14.1 reduct ion program via 
nonlinear transform ation formulae (iHarmanec et al.l 1 19941 
IHarmanec & Homiri 998). 

5. AS AS V magnitude observations were extracted from the 
public AS AS database (Poimanski 2002); we used the data 
from the diaphragm, which gives the smallest rms errors and 
omitted a few clearly deviating data points. 

The journal of all photometric observations is in Table lA.ll 
Homogenized UBV magnitudes of all comparison and check 
stars used can be found in Table lA.2l 

Appendix B: Spectroscopy 

The journal of all spectroscopic observations can be found in 
Table lB.l1 The individual data files are identified there by letters. 
Below, we provide a few comments on them. 

- File A CCD spectra of V2368 Oph covering the wavelength 
region 6260-6750 A. They were secured with a SITe-005 
800 X 2000 CCD detector attached to the medium 0.7-m 
camera of the coude focus of the Ondfejov 2.0 m telescope 
(OND). The spectra were obtained between June 2007 and 
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Table A.l. Journal of available photometry of V2368 Oph 



Station 


Time interval 


No. of 


Passbands 


HD of comparison 


Source 




(HJD-2400000) 


obs. 




/ check star 




1 


54273.4-54977.6 


423 


UBV 


154660/154895 


this paper 


30 


52065.8-54276.8 


152 


UBV 


*) 


this paper 


66 


52765.4-52766.6 


10 


UBV 


154660/154895, 154445 


this paper 


61 


47912.6^9061.9 


111 


V 


all-sky 


Perrvman&ESAd997) 


93 


53055.9-53290.5 


47 


V 


all-sky 


Poimanski ('2002) 



Individual observing stations are distinguished by running numbers they have in the Prague/Zagreb photometric archives — see column "Station"; 01 . . . Hvar 0.65-m, Cassegrain 
reflector, EM19789QB tube; 30 . . . San Pedro Maitir, 0.84-m reflector, Cuenta-pulsos photometer; 61 . . . Hipparcos all-sky Hp photometry transformed to Johnson 1^; 66 . . . TNO 0.40-m 
Cassegrain reflector, SSP5A photometer; 93 . . . ASAS data archive itPoimanskik200Z^ . 

*) All-sky photometry or differential photometry relative to various comparisons during the first season when V2368 Oph was used as a comparison for obseiTations of U Oph, then 
relative to HD 154660 — see the text for details. 



June 2008 and have a linear dispersion of 17.2 A mm"'(red) 
and a 2-pixel resolving power of about 12600 (11-12 km s"' 
per pixel). Their S /N ranges from 50 (1 spectrum) to 370, 
and the majority have S /N over 200. 

- File B CCD spectra covering the wavelength region 6150- 
6750 A with a resolution of 6.6 km s ' per pixel. They were 
obtained at the DAO 1 .22-m telescope between August 2007 
and September 2009 and have a reciprocal linear dispersion 
of 10 A mm^'. The detector used was a SITe-4 4096 x 2048 
CCD, and the 3-pixel resolving power was about 15000. 
Their S/N ranges from 100 to 370. 

- File C CCD echelle spectra secured with the 2. 14-m reflector 
of the SPM obsei-vatory in April 2003. The CCD detector has 
1024 X 1024 pixels, and the setting used covered the wave- 
length region from 4000 to 6700 A in grating orders 33 to 
60. The nominal resolution of the spectrograph is 18000 at 
5000 A, which translates into 2-pixel resolution of about 17 
km s"'. 

- File D Another set of echelle CCD spectra from SPM, se- 
cured in April 2007 

- File E The third set of echelle CCD spectra from SPM, se- 
cured from May 30 to June 1, 2007. The S/N of the SPM 
spectra ranges between 120 and 500 for the red, and 85 to 
230 for the blue parts of the spectra. 

The initial reductions of the DAO spectra (bias subtraction, 
flatfielding and conversion to 1-D images) were carried out by 
SY in IRAF. The initial reduction of the SPM and OND spectra 
was carried out by PE and by Dr. M. Slechta, respectively, in- 
cluding the wavelength calibration. The remaining reductions of 
all spectra (including wavelength calibration for the DAO spec- 
tra, continuum rectification, and removal of c osmics and flaws) 
was carried out with the program SPEFO (iHom et al.l Il996t 
[Skodal[T996l) . 



